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where  T  =4'm\L/[\n(L/2R)  +  0.84]  is  the 
Stokes  drag  coefficient  for  a  cylinder  (26), 
and  Tj  is  the  effective  viscosity  of  the  nematic. 
The  solid  lines  in  Fig.  4B  are  the  results  of 
fits  to  the  solution  of  this  equation 

/  8 Mg  \ 

z{t)  =  -4vKC(Mdzy)  [1  -  exp(-t/T)] 

(6) 

with  the  time  constant,  t  =  TI\4itKC(d^l 
dz)2],  constrained  to  have  the  same  value,  t  = 
270  s,  in  all  three  fits.  These  results  give  clear 
qualitative  support  for  Eqs.  4  to  6  in  describ¬ 
ing  the  levitation.  Quantitative  comparisons 
are  made  difficult  by  the  complicated  aniso¬ 
tropy  of  the  nematic’s  viscoelasticity  (21). 
However,  selecting  characteristic  values  for 
5CB  of  ti  ~  1  Poise  (27)  and^~  5  X  10“7 
dynes  leads  to  t  ^  140  s,  which  is  in  reason¬ 
able  agreement  with  the  measured  value. 

This  manipulation  of  wires  with  director 
gradients  suggests  methods  for  self-assem¬ 
bly  using  more  elaborate  spatially  depen¬ 
dent  director  fields.  For  example,  as  a 
straightforward  extension  of  the  levitation, 
a  liquid  crystal  that  twists  through  many 
periods  (such  as  a  cholesteric  liquid  crystal) 
could  serve  as  host  to  the  creation  of  a 
periodic  array  of  planes  of  oriented  nano- 


Carbon  nanotubes  have  been  considered  lead¬ 
ing  candidates  for  nanoscale  electronic  applica¬ 
tions  (1,  2).  Previous  measurements  of  nano¬ 
tube  electronics  have  shown  electron  confine¬ 
ment  (quantum  dot)  effects  such  as  single¬ 
electron  charging  and  energy-level  quantization 
(3-5).  Nanotube  properties  such  as  long  spin 
lifetimes  make  them  ideal  candidates  for  spin- 
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wires.  Another  possible  application  of  these 
forces  involves  wire  sorting  through  the 
repulsion  from  a  flat  surface.  As  Eq.  3 
indicates,  gravity  and  the  repulsive  force 
balance  at  a  height  that  varies  to  first  order 
with  R~2  and  is  independent  of  L  when  0  is 
fixed;  for  example,  by  a  large  external 
field.  Hence  this  repulsion  can  be  used  to 
spatially  separate  wires  of  different  radii. 
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The  device  under  study  consists  of  a  nano¬ 
tube  —2  nm  in  diameter,  of  length  L  —1.5  |xm 
between  metal  contacts,  with  three  top  gates 
(11, 12)  and  a  doped  Si  backgate  (Fig.  1)  (13). 
Room  temperature  measurements  of  conduc¬ 
tance  as  a  function  of  any  gate  voltage  indicate 
either  that  the  nanotube  is  metallic  with  a  strong 
resonance  or  that  it  is  a  small  band-gap  semi¬ 
conductor.  At  low  temperatures,  single-electron 
charging  is  observed  as  Coulomb  blockade  con¬ 
ductance  peaks,  measured  as  a  function  of  ap¬ 
plied  gate  voltages.  The  observed  peak  pattern 
is  different  for  each  gate  and  in  some  regions 
shows  considerable  superstructure,  which  ap¬ 
pears  with  increasing  source-drain  voltage,  VSD. 
We  observe  a  double  periodicity  of  Coulomb 
charging  phenomena,  which  is  consistent  with  a 
nanotube  quantum  dot,  defined  by  tunnel  bar¬ 
riers  to  the  leads,  that  has  been  split  into  two 
dots  of  roughly  equal  size  by  a  defect  in  the 
middle  of  the  device  (14).  Structural  defects 
often  create  tunable  tunnel  barriers  (15-18)  and 
typically  appear  in  nanotubes  longer  than  —200 
nm.  Although  such  defects  can  be  controllably 
fabricated  (19),  in  this  case  the  defect  is  inher¬ 
ent  to  the  tube.  The  resulting  dots  each  have  a 
quantized  energy  level  spacing  of  A  =  hvF/ 
2L~l  meV  (assuming  only  spin  degeneracy), 
where  vF  =  8.1  X  105  m/s  is  the  Fermi  velocity 
for  a  metallic  nanotube  (4). 

The  top  gate  voltages  have  differential  ca¬ 
pacitive  coupling  to  the  quantized  energy  levels 
of  the  two  dots.  Transport  measurements  indi- 
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We  have  measured  carbon  nanotube  quantum  dots  with  multiple  electro¬ 
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interactions  in  the  quantum  regime  in  a  molecular  conductor  is  important 
for  applications  such  as  quantum  computation. 
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cate  that  gate  G1  predominantly  couples  to  dot 
1  and  gate  G2  to  dot  2.  In  the  transconductance, 
dI/dVGl,  as  a  function  of  VG1  and  VG2  (Fig. 
2A),  paired  regions  of  large  dI/dVGl  appear  on 
the  vertices  of  a  hexagonal  lattice  (13).  The 
resulting  honeycomb-shaped  array  of  high  dl/ 
dVG1  can  be  contrasted  with  the  straight  lines 
that  would  appear  if  each  applied  gate  voltage 
coupled  equally  to  the  dots.  The  observed  pat¬ 
tern  strongly  resembles  the  charge  stability  di¬ 
agram  found  for  weakly  tunneling  double  quan¬ 
tum  dot  (interdot  tunneling  <$C  2e2/h),  where 
each  cell  of  the  honeycomb  corresponds  to  a 
well-defined  electron  configuration  for  the  dou¬ 
ble  dot  (20).  The  vertex  pairs  are  points  of 
degeneracy  between  the  two  dots  that  have  split 
because  of  interdot  coupling  (either  capacitive 
or  tunneling).  At  these  vertices,  energy  levels 
resonantly  align  between  the  dots  as  well  as  the 
leads  (Fig.  2C).  For  a  series  double  dot,  reso¬ 
nant  transport  only  occurs  at  the  vertices;  finite 
nonresonant  conductance  along  the  edges  of  the 
honeycomb  cells  (Fig.  2,  D  and  E)  occurs  via 
cotunneling,  which  is  suppressed  in  the  weak- 
tunneling  regime  (20). 

At  finite  bias,  the  vertex  points  expand  into 
triangular  regions  (Fig.  2B)  defined  by  lines 
where  the  energy  levels  of  each  dot  align  with 
either  the  top  or  bottom  of  the  bias  window. 
According  to  theory,  resonant  tunneling  should 
occur  only  along  the  one  side  of  the  triangle  that 
connects  the  original  vertices.  The  observation 
of  off-resonance  tunneling  in  the  body  of  the 
triangles  has  previously  been  attributed  to  both 
inelastic  phonon  contributions  and  cotunneling 
processes  (21-23).  Because  it  appears  that  co¬ 
tunneling  is  suppressed  in  this  region  (we  do 
not  see  conduction  along  the  edges  of  the  hon¬ 
eycomb),  conductance  inside  the  finite-bias  tri¬ 
angles  suggests,  perhaps,  a  substantial  contribu¬ 
tion  from  electron-phonon  coupling  to  trans- 
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Fig.  1.  (A)  Schematic  of  top-gated  device.  (B) 
Electron  micrograph  of  a  representative  device. 
Arrows  indicate  the  embedded  nanotube. 


port.  However,  this  conclusion  appears  incon¬ 
sistent  with  the  predicted  suppression  of 
electron-phonon  coupling  in  nanotubes  at  low 
temperatures  ( 24 ,  25).  It  seems  more  likely  that 
inelastic  processes  are  occurring  through  elec¬ 
tron  interactions  with  the  source/drain  reser¬ 
voirs  or  the  substrate.  Figure  2B  does  not  show 
resonant  excited-state  levels  inside  the  finite- 
bias  triangles,  as  are  observed  in  lateral  quan¬ 
tum  dots  (20).  This  is  not  surprising,  as  eVSD  < 
A  in  the  present  case.  Increasing  VSD  causes  the 
triangles  to  grow,  as  expected,  but  also  inhibits 
our  ability  to  resolve  individual  vertices  and 
even  individual  honeycomb  cells.  Overall, 
however,  the  transport  characteristics  of  this 
nanotube  double  dot  are  qualitatively  similar  to 
those  observed  for  double  dots  fabricated  in 
semiconductor  heterostructures  (27,  26).  Such 
similarities  suggest  that  more  complex  quantum 
devices  formed  from  molecular  conductors 
such  as  nanotubes  will  continue  to  share  fea¬ 
tures  with  well-studied  semiconductor  quantum 
dot  systems.  However,  the  larger  energy  scales 
in  carbon  nanotubes  allow  operation  at  higher 
temperatures. 

Double  dot  capacitances  are  obtained  by 
measuring  the  sizes  and  separations  of  the 
honeycombs  and  vertices  (Fig.  2)  (20).  The 
size  of  the  honeycomb  in  Fig.  2A  determines 

AFG1(2)>  Where 


AF01(2)  =  ?i^(l  +  A  /Ecm)  (1) 

^Gl(2) 

with  CG1(2)  the  capacitance  of  dot  1(2)  to 
gate  1(2)  and  £C1(2)  ~  2meV  the  charging 
energy  for  dot  1(2).  Ec i(2)  is  determined 
from  \qVsd\  measured  at  the  apex  of  a  con¬ 
ductance  diamond  in  a  plot  of  conductance 
as  a  function  of  VSD  and  any  gate  voltage 
(not  shown).  From  the  size  of  the  finite-bias 
triangles  in  Fig.  2B,  we  determine  capaci¬ 
tance  ratios  using 

C  \(2) 

§VGl(2)  =  7^1VSD  (2) 

UG  1(2) 

With  the  above  relationships,  we  calculate 
typical  dot  capacitances  of  C1(2)  ~  160  aF 
and  CG1(2)  ~  45  aF. 

The  amount  of  capacitive  coupling  and 
interdot  tunneling  can  be  determined  from 
the  vertex  splitting  (27,  28).  These  cou¬ 
plings  cause  an  energy  shift  related  to  the 
formation  of  charge  polarized,  or  “bond¬ 
ing,”  states  between  the  quantum  dots 
(tunnel  coupling  is  analogous  to  covalent 
bonding,  and  capacitive  coupling  to  ionic 
bonding).  Both  coupling  parameters  can  be 
tuned  with  gate  voltages  that  open  the  in¬ 
terdot  constriction,  so  it  is  difficult  to 
determine  their  relative  contributions  to 


Fig.  2.  Experimental  charge  sta¬ 
bility  diagrams  for  the  series 
double  quantum  dot  as  a  func¬ 
tion  of  two  gate  voltages,  each 
shifting  the  energy  levels  of  a 
single  dot.  Voltages  VC1  and  VG2 
are  divided  by  five  before  being 
applied  to  the  gates  of  the  de¬ 
vice.  (A)  Color  scale  displays  dl/ 
dVc1  calculated  from  dc  current 
(lDC)  at  Vs^  =  500  [xV.  White 
lines  are  guides  to  the  eye  show¬ 
ing  the  honeycomb  pattern  of 
peaks  in  conductance.  Vertex 
pairs  correspond  to  points  of  de¬ 
generacy  between  the  two  dots 
where  resonant  transport  occurs, 
whereas  cotunneling  may  pro¬ 
duce  finite  conductance  along 
the  honeycomb  edges.  (B) 

Zoom-in  of  a  vertex  pair  at  VSD  = 

500  |jlV  with  lines  of  constant 
energy  designated  by  white 
dashed  lines.  Colorscale  dis¬ 
plays  lDC.  Vertex  dimensions 
are  indicated  on  the  diagram, 
as  are  particular  electron  con¬ 
figurations  on  the  two  dots  for 
different  regions  of  gate  space. 

Note  the  triangular  shape  re¬ 
sulting  from  finite  VSD,  and 
nonzero  off-resonance  conduc¬ 
tance  within  the  triangles.  (C) 

Level  diagrams  for  a  double 
quantum  dot  depicting  reso¬ 
nant  transport  with  Tm  the  in¬ 
terdot  tunnel  rate.  The  bias  window  is  VSD=  ([Ls-\xD)/e.  (D  and  E)  Level  diagrams  depicting 
configurations  where  cotunneling  may  contribute  to  transport. 
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vertex  splitting.  If  we  assume  the  vertex 
splitting  is  due  to  capacitive  coupling 
alone,  we  find  a  mutual  capacitance  be¬ 
tween  dots  Cm  ~  90  aF,  using 


A*Sip)  = 


\e\Cm 

Cg1(2)C2(1) 


(3) 


where  is  the  horizontal  or  vertical 

component  of  the  diagonal  splitting  measured 
between  vertices  (20).  A  separate  determina¬ 
tion  of  the  mutual  capacitance  based  on 
single-dot  capacitances  (determined  above) 
and  the  magnitude  of  the  charging  energy 
using  the  relation 

e 2 

ECH2)  =  —  (1  -  C^/C^r1  (4) 

'“'1(2) 


gives  Cm  —100  aF.  The  similarity  between 
these  two  values  for  Cm  indicates  that  the 
vertex  splitting  in  this  region  is  dominated  by 
capacitive  coupling.  The  interaction  energy  is 
then  approximately 

e1 

Em  =  —  (C,C2/C2m-  l)-1  (5) 

^ m 

-  700  |xeV. 

Opening  the  interdot  constriction  is  ex¬ 
pected  to  increase  the  tunnel  coupling  expo¬ 
nentially  faster  than  the  capacitive  coupling 
(27).  We  can  therefore  compare  the  amount 
of  interdot  tunnel  coupling  in  different  gate 
configurations  using  the  fractional  peak  split¬ 
ting,  defined  f=2AVs/AVP  where  AVS  is  the 
diagonal  splitting  measured  between  vertices 


Fig.  3.  (A)  Experimental  charge 
stability  diagram  demonstrating 
weak  interdot  tunnel  coupling. 
VSD  =  500  |mV.  Vertex  peak  split¬ 
ting  (AVS)  and  honeycomb  cell 
dimensions  (AVp)  indicated  on 
the  plot.  (B)  Moving  to  a  differ¬ 
ent  region  of  VC1  changes  the 
tunnel  coupling  between  dots, 
producing  a  smeared  charge  sta¬ 
bility  diagram.  (C  to  E)  Experi¬ 
mental  stability  diagrams  for 
fixed  charge  configuration  at 
various  VBG.  VSD  =  800  pV.  The 
capacitive  coupling  of  the  back- 
gate  to  the  two  dots  requires 
that  slightly  different  ranges  of 
VC1  and  VC2  be  used  to  cover  the 
same  charge  configurations  at 
different  VBG  values. 


Fig.  4.  Temperature 
dependence  of  peak 
heights  for  two  differ¬ 
ent  regimes  of  trans¬ 
port.  One-dimensional 
sweeps  taken  in  VG2 
with  VC1  =  -1800  mV. 

Representative  cotun- 
neling  peak  (A),  de¬ 
noted  "Cot."  and  ver¬ 
tex  peaks  (B),  denoted 
"Res.,"  at  tempera¬ 
tures  between  T 
350  mK  and  915  mK. 

(C)  Peak  height  as  a 
function  of  tempera¬ 
ture  on  a  log-log  plot. 

Peaks  are  sorted  by 
type  (resonant  or  co¬ 
tunneling)  into  two 
distinct  groups  and 
have  all  been  offset  to 
collapse  approximate¬ 
ly  onto  a  single  curve 
for  each  group.  The 
black  line  shows  the 
average  power  law  be¬ 
havior,  T1-3,  for  cotunneling  peaks.  Peaks  shown  in  (A)  and  (B)  are  shown  in  red.  (Inset) 
Representative  honeycomb  cell  in  the  region  of  measurement. 
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and  AVP  is  the  distance  between  vertex  pairs 
(Fig.  3A)  (26).  A  value  off  =  1  corresponds 
to  a  transparent  interdot  constriction,  while 
/  =  0  indicates  two  isolated  dots  in  the  limit 
of  zero  mutual  capacitance  between  dots. 
Changes  in  f  however,  should  be  dominated 
by  tunnel  coupling,  so  measuring  how  / 
changes  as  a  function  of  gate  voltage  allows 
the  contribution  to  the  vertex  splitting  from 
interdot  tunneling  to  be  inferred.  We  observe 
an  evolution  of /by  holding  the  range  of  VG2 
fixed  and  examining  different  regimes  of  VGl 
(Fig.  3,  A  and  B).  In  Fig.  3  A,  where  honey¬ 
comb  cells  and  vertices  are  well  defined,  we 
find  that  /  —  0.3.  This  regime  can  be  con¬ 
trasted  with  that  of  Fig.  3B,  where  an  increase 
in  interdot  tunnel  coupling  leads  to  a  smear¬ 
ing  of  vertices  and  honeycomb  cells.  In  this 
region,  vertex  splitting  grows  to  a  value  of 

f~  0.7. 

In  Fig.  3,  C  to  E,  we  use  a  combination  of 
VGl  and  VG2  to  maintain  a  fixed  electron  con¬ 
figuration  while  tuning  tunneling  with  the  back- 
gate.  The  evolution  evident  here  shows  how 
varying  VBG  from  250  mV  to  375  mV  changes 
our  stability  diagram  from  weak  to  strong  co¬ 
tunneling  regimes.  The  corresponding  increase 
in  vertex  heights  indicates  that  we  are  also 
influencing  the  total  amount  of  tunneling 
through  the  double  dot.  Because  the  vertex 
splitting  does  not  change  substantially,  it  seems 
that  the  backgate  is  primarily  tuning  the  trans¬ 
parency  of  the  tunnel  barriers  at  the  source 
and/or  drain  electrodes  (15). 

Increasing  temperature  increases  cotunnel¬ 
ing  conductance  (i.e.,  conductance  on  honey¬ 
comb  edges)  by  a  relation  well  described  by  a 
power  law,  Gcot<*Tr[,  over  the  range  investigat¬ 
ed,  where  t|  is  between  1  and  2  with  a  typical 
value  ofri~1.3  (Fig.  4,  A  and  C,  upper  curve). 
In  contrast,  resonant  vertex  peak  conductance  is 
roughly  independent  of  temperature,  or  only 
slightly  increases  with  increasing  temperature 
(Fig.  4,  B  and  C,  lower  curve).  Power  law 
behavior  with  r\  =  2  has  been  predicted  for 
cotunneling  peaks  in  the  weak  interdot  tunnel¬ 
ing  regime  (Gint  <$C  e2/ h)  at  low  temperatures 
(kBT  <$C  A)  (29).  A  power  law  with  r\  = 
1.25 — comparable  to  what  we  find — is  predict¬ 
ed  (29)  for  intermediate  tunnel  coupling,  and  A 
<$C  kBT  (corresponding  to  temperatures  higher 
than  those  of  our  measurements).  The  nearly 
temperature-independent  behavior  we  observe 
for  resonant  peaks  does  not  appear  consistent 
with  predictions  of  peak  height  scaling  as  1/T  in 
the  low-temperature  regime  (28, 30).  Given  that 
other  measurements  of  this  system  correspond 
relatively  well  to  predicted  behavior,  it  is  sur¬ 
prising  that  the  temperature  dependence  devi¬ 
ates  so  far  from  theoretical  expectations. 
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recent  ab  initio  study  of  the  aqueous  liquid- 
vapor  interface  using  32  waters  produced  inter¬ 
esting  results  regarding  dangling  bonds  but 


Fig.  1.  Snapshot  of  the  calculated  aqueous  liq¬ 
uid-vapor  interface.  The  individual  water 
monomers  are  represented  by  yellow  cylinders. 
The  colored  isosurfaces  on  the  top  represent 
the  HOMO,  whereas  the  colored  isosurfaces  at 
the  bottom  represent  the  LUMO.  The  orbitals 
are  obtained  by  direct  diagonalization  of  the  KS 
matrix,  indicating  that  the  reactivity  of  the 
water  slab  is  localized  on  the  surface  for  pro¬ 
tons  and  electrons. 


An  ab  Initio  Molecular  Dynamics 
Study  of  the  Aqueous 
Liquid-Vapor  Interface 

l-Feng  W.  Kuo  and  Christopher  J.  Mundy* * 

We  present  an  ab  initio  molecular  dynamics  simulation  of  the  aqueous  liquid- 
vapor  interface.  Having  successfully  stabilized  a  region  of  bulk  water  in  the 
center  of  a  water  slab,  we  were  able  to  reproduce  and  further  quantify  the 
experimentally  observed  abundance  of  surface  “acceptor-only”  (19%)  and 
“single-donor”  (66%)  moieties  as  well  as  substantial  surface  relaxation  ap¬ 
proaching  the  liquid-vapor  interface.  Examination  of  the  orientational  dynamics 
points  to  a  faster  relaxation  in  the  interfacial  region.  Furthermore,  the  average 
value  of  the  dipole  decreases  and  the  average  value  of  the  highest  occupied 
molecular  orbital  for  each  water  molecule  increases  approaching  the  liquid- 
vapor  interface.  Our  results  support  the  idea  that  the  surface  contains,  on 
average,  far  more  reactive  states  than  the  bulk. 
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